Mon. Not. R. Astron. Soc. 000,[T]|^(2004) Printed 2 February 2008 (MN MbX style file v2.2) 



Reflection spectra from an accretion disc ifluminated by a neutron 
star X-ray burst 



D. R. Ballantyne^ 



O 
O 

(N 



> 

o 

\o 

(N 
O 

o 
o 



Canadian Institute for Theoretical Astrophysics, 60 St. George Street, Toronto, Ontario, Canada M5S 3H8 



2 February 2008 



ABSTRACT 

Recent time-resolved X-ray spectra of a neutron star undergoing a superburst revealed an 
Fe Ka line and edge consistent with reprocessing from the surrounding accretion disc. Here, 
we present models of X-ray reflection from a constant density slab illuminated by a blackbody, 
the spectrum emitted by a neutron star burst. The calculations predict a prominent Fe Ka line 
and a rich soft X-ray line spectrum which is superimposed on a strong free-free continuum. 
The lines slowly vanish as the ionization parameter of the slab is increased, but the free-free 
continuum remains dominant at energies less than 1 keV. The reflection spectrum has a quasi- 
blackbody shape only at energies greater than 3 ke V. If the incident blackbody is added to the 
reflection spectrum, the Fe Ka equivalent width varies between 100 and 300 eV depending 
on the ionization parameter and the temperature, kT, of the blackbody. The equivalent width 
is correlated with kT, and therefore we predict a strong Fe Ka line when an X-ray burst is 
at its brightest (if iron is not too ionized or the reflection amplitude too small). Extending 
the study of reflection features in the spectra of superbursts to lower energies would provide 
further constraints on the accretion flow. If the Fe Ka line or other features are relativistically 
broadened then they can determine the system inclination angle (which leads to the neutron 
star mass), and, if the mass is known, a lower-limit to the mass/radius ratio of the star 
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1 INTRODUCTION 

The reprocessing of external X-rays by accreting material has 
been observed from active galactic nuclei (AGN) and Galac- 
tic black hole candidates (GBHCs) for over a decade (e.g., 
[founds et al. 199 0; Nandra & Pounds 1 9 94 ;_ Gierlinski et al,,, 
ri999l: fBallantvne. Vaughan & Fabiarjj l200i iMiller et all boolT 



In these objects, the origin and location of the illuminating 
X-ray power-law is largely unknown, although it is widely 
believed to be either within a m a gnetic ac cretion disc corona 
(e.g., |Gal£eV j_^osngr Vaianal |l97^ iHaardt &_MMaschi 
1993 ; iHaard t. Maraschi & Ghisellini 1994^ iDi Mattec 
1998 ) or a centrally located geometrically thick flow (e.g., 



Sha^iro^Lightman & Eardlev 1976; Zdziarski, Lubinski & Smith 



1999; Zdziarski et al. 2003). Nevertheless, many models of X-ray 
reflection from black hole accretion discs have been computed 



jGeorge & Fabianll99lLlMatt Perola & Piroll99lHRoss & Fabianl 
19931 IZ vckietal.' '1994*: 'Ross, Fabian & Young' '1999'; 



Navakshi n. Kazanas & Kallman 2000; Navakshin & Kallman 
2001; Ballantvne, Ross & Fabian 200 Ibt iRozanska et alJ l200l 
Ballantvne. Turner & Blaes. 2004), and, in some cases, s uccess- 
fully applied to data teallantvne. Iwasawa & Fabiai] 12001 at 
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l Orr et alJ l200lt iDe R osa et alj l200i iLonginotti et alj l2003t 

Barri o. Done & Navak shin 2003: lMiller et alj2004 ! 

Recently, IStrohmaver & Brownl i2002l) discovered an Fe Ka 
line and edge in the X-ray spectra of the low-mass X-ray binary 
(LMXB) 4U 1820-30 as it was undergoing a superburst — a pow- 
erful thermonuclear explosion on the surface of the neutron star that 
can la s t many hours fc ornelisse et al. 2000; Strohinaver & Bildst^ 
l2003t iKuulkersI 120031 Based on the strength of the Fe Ka 
line, these authors suggested that the features were caused by 
reflection of the burst bl ackbody spec t rum f rom the surround- 
ing accretion disc (e.g., iDav & Dond Il99ll) . This hypothesis 
was supported by spectral fitting with detailed reflection models 
(Ballantvne & Strohmavei 2004). The unique aspect of reflection 
during a burst from a LMXB is that there is no uncertainty in the lo- 
cation of the X-ray source, as it is located on the surface of the neu- 
tron star and outshines the persistent X-ray emission from the disc. 
This fact allows for far less ambiguity when analyzing changes to 
the reflection features , as they are more likely to trace the evolution 
in the accretion flow iBallantvne & StrohmaveJ2004h . 

X-ray emission lines have also been observed in the persis- 
tent emission f r om LMXBs (e.g..[S male et alll993l : lAngelini et alJ 
119951: ISchulJ Il999t lAsai et alj l2000h . and these are also 
thought t o be caused by reprocess i ng by matter in the accre- 
tion flow iShakura & SunvaeMll973t iMever & Mever-HofmeisteJ 
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Il982h . The spectral shape of the persistent emission is con- 
sistent with bremsstrahlung with a tem perature of 5-10 keV 
iLiu. van Paradiis & van den Heuvell 12001 'I . most likely arising 
from the boundary layer between the disc and the neutron star De- 
tailed models of the photoionized layer on the disc have been per- 
formed and predict a wealth of emission lines in the soft X-ray band 
(e.g ., Ko & Kallmari ll994l : Ijimenez-Garate. Raymond & Liedahll 
[2002) that can be compared against observations made with the 
grating spectrometers onboard Chandra and XMM-Newton ( e.g., 
ISchulz etalJboOll ICottam et alJl200lt iKallman et all l2003h. In 
contrast, a X-ray burst e mits a black body spectrunT TswMik et alJ 
[l977; Hoffman, Lewin & Dotvlll977l) . and can reach the Edding- 
ton luminosity (LEdd ~ 10''* erg s~^), well over an order of 
magnitude greater than the pe rsistent luminosity of many LMXBs 
iStrohmaver & Bildstenl2003l) . Thus, disc reflection will be domi- 
nated by the blackbody component during the burst. 

This paper presents reflection spectra from a uniform accre- 
tion disc illuminated by a blackbody, as in a LMXB during a 
Type I X-ray (super)burst. T hese spectra were used to fit the super- 
burst data of 4U 1820-30 jBallantvne & Strohmave3l2004h . but, 
as shown below, they contain much more information than what 
was actually needed for those Rossi X-ray Timing Explorer (RXTE) 
data. Observations of reflection features from Type I bursts or su- 
perbursts with sensitive, large bandwidth instruments such as the 
X-ray Telescope (XRT) on Swift could provide a wealth of new in- 
formation on the structure and behaviour of accretion discs. 

In the next section we outline the reflection calculations, and 
present the resulting spectra in Section |3| We conclude by dis- 
cussing the results in Section|4] 



2 CALCULATIONS 

The reflection spectra were computed using the code of 
IRoss & Fabian ( 1993.) . Details of the computational tec hniques can 
be fou nd in that paper and in referen ces therein (see also lRoss et alj 
^1999^) and'Ba llantvne et alj|200ld) . and so only a brief overview 
will be presented here. 

A one-dimensional slab of material with a constant hydrogen 
number density nn is illuminated by a blackbody spectrum with 
temperature kT and flux (defined between 1 eV and 100 keV) Fx. 
Thus, we define the ionization parameter 



The input radiation is tra nsferred into the layer us ing the 
two-stream approximation tevbicki & Lightman 'l979j) while 
the outgoing radiation is treated using the modified Fokker - 
Planck/diffusion method of iRoss. We aver & McCra^ h978l) . 
Therefore, the output spectrum is averaged over all viewing an- 
gles. The ionization and thermal structure of the illuminated gas 
is iterated until each zone reaches ionization and thermal bal- 
ance, and at this point the final reflection spectrum is computed. 
The illuminated slab always has a large enough Thomson depth 
(tt =6-20, depending on 5) so that the highest energy photons 
wi ll interact many times w i th the gas. The cosmic abundances 
of|M omson & McCammor] h983h were assumed, but hydrogen 
and helium are taken to be completely ionized at all times (i.e., 
rie = 1.2nH). The following ionization stages of the most impor- 
tant astrophysical metals are treated in the calculation: C V-VII, 
N VI-VIII, O V-IX, Mg IX-XIII, Si XI-XV and Fe XVI-XXVII. 

Models were calculated for 1 ^ kT ^3.15 keV, as this spans 



the range most often observed during X-ray bursts jKuulkersl2003l) . 

and for 1 ^ log ^ ^ 3.75. To change 5, tih was fixed at 10^^ cm~^, 
and Fx was varied. 



3 RESULTS 

We begin by considering t he evolution of t he reflection spectra as a 
function of ^ (cf. Fig. 2 in R^ssetai]l99i). FieurelTlshows a series 
of spectra for four different values of kT. Each panel shows reflec- 
tion spectra with log ^ = 1,1.5,2, 2.5, 3, 3.5 (from bottom to top), 
thereby spanning the range from very weakly ionized reflectors to 
highly ionized ones. For the log ^ = 2 model, we have also plotted 
the incident blackbody spectrum. 

Unlike the power-law continua of AGN and GBHCs where the 
radiant energy per frequency decade reaches a maximum at ener- 
gies > 10 keV (for photon-indices F < 2) where atomic absorption 
is minimal, the blackbody spectra peak at energies where the pho- 
toelectric absorption by metals is very important. As a result, only 
a small fraction of the photons incident on the slab are scattered 
back without first being absorbed by an ion. The reflection spectra 
therefore only mimics a blackbody at energies E > 3 keV, with 
the agreement increasing as log ^ grows and the metals are ionized. 
The energy absorbed by the metals in the gas is thermalized and 
reemerges as a strong soft free-free continuum with a variety of re- 
combination lines superimposed. As this thermal emission includes 
the power absorbed at higher energies, it outshines the incident 
blackbody at _E < 1 keV for all values of ^ considered. The rich 
emission line spectrum is reminiscent of those predicted from mod- 
els of reprocessing of the persistent emission ( Jimenez-Garate et al] 
[2002). However, in the most luminous stage of an X-r ay burst, the 
disc may be severely ionized jBallantvne & Strohmaver>i2004.) . and 
will leave a Compton broadened O VIII Lya line as the only emis- 
sion feature in the soft X-ray region of the spectrum. 

It is instructive to look at the most important heating and cool- 
ing processes within the illuminated layer A plot of the heating and 
cooling rates as a function of Thomson depth into the slab for the 
log ^ = 3, kT — 2.5 keV model is shown in the lower panel of Fig- 
ure[2| In this figure, the heating and cooling processes are shown by 
red and blue lines respectively, while the different line styles sepa- 
rate the physical process (for example, Compton heating is denoted 
by the red dotted line, but Compton cooling is shown with the blue 
dotted line). The upper panel of Fig.[2|shows the final iron ioniza- 
tion structure in the slab. Hydrogenic iron dominates between the 
surface and the first Thomson depth, but recombines to He-like Fe 
before reaching tt = 2. The iron recombines rapidly just beyond 
this point and reaches the lowest ionization stage treated (Fe XVI) 
at Tt ~ 5. Turning to the heating and cooling rates, we see that 
free-free emission is the dominant cooling mechanism over the en- 
tire depth of the slab. At the surface, where the gas temperature 
reaches T ~ 5 x 10® K (the black line in the lower panel of Fig.[2j, 
the bremsstrahlung cooling is balanced by Compton and photoelec- 
tric heating. However, the Compton heating rate falls rapidly with 
Tt because there are very few high energy photons in the incident 
blackbody spectra. In contrast, the photoelectric heating rate in- 
creases with Tt until it reaches a maximum at tt ~ 3, where iron 
is rapidly recombining (upper panel of Fig.[2j and the photoioniza- 
tion rates jumps. The line cooling rate also increases rapidly here, 
but cannot overcome the free-free cooling. After the iron has al- 
most fully recombined at tt ~ 5, the photoelectric heating rate 
falls (very few ionizing photons can penetrate to such depths) and 
the dominant heating process becomes free-free absorption. The 
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Figure 1. Examples of reflection spectra from an accretion disc illuminated by a blackbody continuum. The four panels show results for different values of kT. 
Within each panel 6 reflection spectra are shown with solid lines. They denote models with (moving from bottom to top) log § = 1, 1.5, 2, 2.5, 3, 3.5. All of the 
spectra, except for the log 5 = 2 model, have been vertically offset from one another for clarity. The dotted line in each panel shows the illuminating blackbody 
for the log ^ = 2 model. 



line-cooling rate jumps one more time at tt ~ 9 where oxygen is 
recombining (note the corresponding fall in the gas temperature). 
The dominance of free-free cooling and photoelectric heating in 
the thermal structure of the illuminated slab emphasizes the dom- 
inance of the absorption-thermalization process in these models. 
The observational outcome is the rich X-ray line spectrum seen be- 
low 1 keV in FigureQ 

Similar to the reflection spectra predicted for the environment 
of accreting black holes, the Fe Ka line is the most important spec- 



tral feature above ~ 3 keV in the blackbody models, and is an 
important indicator of the ionization parameter of the reflecting 
medium. Fig. Q shows that in each panel the Fe Ka line evolves 
from a strong 6.4 keV line to a weak, Auger-destroyed line, and 
then to a broad 6.7 keV feature before disappearing entirely due to 
the complete ionization of iron. This is the same sequence as that 
presented by Ross et al. ( 1999) for illumination by a power-law. 

The equivalent width (EW) of the Fe Ka line as a function 
of log ^ is shown in Figure|3| The top panel shows the EWs mea- 
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Figure 2. Details from the log § = 3, kT = 2.5 keV reflection calculation 
as a function of Thomson depth into a slab with njj = 10^^ cm"''. The up- 
per panel plots the iron ionization structure. Moving from the surface down 
into the slab, the iron ions encountered are Fe XXVII (long-dashed line), hy- 
drogenic Fe XXVI (soUd line), Fe XXV (short-dashed line), Fe XVII-XXIV 
(dotted lines) and Fe XVI (solid line; the most neutral iron ion treated). The 
lower panel plots the (H)eating (red curves) and (C)ooling (blue curves) 
rates and the gas temperature (right-hand axis; black solid line). The line 
types distinguish the various physical processes: Compton heating/cooling 
(dotted lines), free-free heating/cooling (short-dashed lines), recombination 
heating/cooling (long-dashed lines), photoelectric heating/line cooling (dot- 
short-dashed lines), and three-body interactions (dot-long-dashed lines). 



sured from the reflection spectra themselves, while the lower panel 
includes the diluting effect of adding in the incident blackbody that 
will also be observed along with the reflection spectrum. The mea- 
surements were made for the four values of kT shown in Fig. Q 
but for 75 values of log ^ between 1 and 3.75. There are two peaks 
in the EW evolution: one when log ^ ~ 1 and Fe Kq is at 6.4 keV, 
and the other when log ^ ~ 2.7 and the line is at 6.7 keV. The EW 
of the cold line is much stronger than the ionized one in the case 
of pure reflection, but is switched when the incident blackbody is 
included. This is because at larger values of ^ there is less metal 
absorption in the illuminated slab and the reflection spectrum is not 
as diluted by the blackbody as in the low ^ models. 

The Fe Ka EW is a strong linear function of the tempera- 
ture of the blackbody; for example, when log ^ = 1 and the inci- 
dent blackbody is included, the EW increases from ~ 70 eV (at 
kT = 1.5 keV) to 180 eV (at fcT = 3 keV). This is a result of 
the increase in photons above the iron K edge (at 7.1-9 keV) as 
kT is increased. This has important implications for the study of 
X-ray bursts, because the blackbody emission is at its hottest when 
the burst is brightest. Therefore, unless iron is fully ionized by the 
burst (i.e., log^ > 3.7), the Fe Kq line should be most prominent 
at the point where the burst is most easily observable. However, a 
superburst may cause changes in the accretion geometry close to 
the neutron star which wo uld reduce the reflection amplit ude and 
the EW of the Fe Ka line jBallantvne & StrohmaveJ2004) . 



Figure 3. The Fe Ka equivalent width (EW) as a function of log § for 
the four values of kT shown in Figure □ The upper panel plots the EW 
measured from the pure reflection spectrum (i.e., there is no contribution 
from the blackbody). The lower panel plots the EW measured from the sum 
of the reflected and incident spectra. In both panels the different line types 
separate the values of kT: 1.5 keV (solid hne), 2 keV (dotted hne), 2.5 keV 
(short-dashed line), and 3 keV (long-dashed line). 

4 DISCUSSION 

Searches for spectral features during X-ray bursts from neutron 
stars have been ongoing for many years. Since the explosion is oc- 
curring on the surface of the star, the spectral features would be 
expected to be redshifted and would therefore be a measure of the 
mass-to-radius ratio of the star (via 1 + 2: = (1 — GM / (? R)~^^^), 
a number needed to constrain the many possible equations of state 
of nuclear matter. During the 1980s there were a few reports of 
X-ray absorption lines at 4.1 keV during Type I X-ray bursts 
dWa ki et alJl984lNakamura. Inoue & TanakJl988l;lMagnier et alJ 
1989i), which were interpreted as redshifted Lya absorption from 
He-like iron. However, the equivalent width of these absorption 
lines were on the order of hundreds of eV, which theoretical 
models of spectral formation in burs t atmospheres could not re - 
produce froster. Ross &. Fabian|[l987l; lOav. Fabian & Rosslll992h . 
Observations by more sensitive instruments on RXTE and Bep- 
poSAX have been unable to discover any other examples of a 
4.1 keV absorption line. Recently, the Reflection Grating Spec- 
trometer (RGS) on XMM-Newton has uncovered evidence for red- 
shifted absorption lines during the X-ray bursts of EXO 0748-676 
( Cottam. Paerels & Mendez 2002), although data from 28 bursts 
were needed in order to obtain the minimum signal-to-noise nec- 
essary to find the lines. 

As first suggested by IPav & Pond il99ll) . X-rays from an 
explosion on a neutron star can be reprocessed by the surround- 
ing accretion disc, leading to features such as an Fe Kq line 
and edge. The analysis of these features can lead to many new 
insights about the structure and evolution of the accretion disc 
( Ballantvne & Strohmaver 2004), as well as basic information 
about the neutron star itself. For example, if the emission lines are 
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relativistically broadened then modeling can re veal the incl ination 
angle of the system to the line of sig ht (FabianetaOElii), a pa- 
rameter often needed to obtain masses of the neutron star and its 
binary companion. Modeling of relativistic lines also produce a ra- 
dius where the reflection originates (assuming some emissivity pro- 
file). This radius is measured in gravitational units {GM/c?), so if 
the mass of the neutron star is known, the radius can be converted to 
physical units, allowing a determination of a mass-to-radius ratio. 
Since the reflection may occur very clos e to the surface of the neu- 
tron s tar (10-20 GM jt? in 4U 1820-30: lBailantvne & StrohmaveJ 
|2004) . the measurement may provide a useful lower-limit to the 
MjR ratio for the neutron star. 



RXTE has so far been the most successful telescope to obtain 
time-resolved spectroscopy of X-ray bursts, especially o f the more 
energetic superbursts (e.g.. lStrohmaver & Brownll2002h . The tele- 
scope produces data for energies greater than 3 keV which allows 
a good determination of the blackbody spectrum and, if detected, 
the Fe Ka line. Yet, the reflection models predict a wealth of infor- 
mation at energies less than 3 keV (Fig.0 that can provide further 
constraints on the abundances, ionization state and dynamics of the 
accretion disc. These features will also strongly depend on the den- 
sity of the disc because of the importance of collisional effects and 
free-free emission/absorption. The spectra shown in Fig. Q were 
calculated assuming a density of tih = 10^^ cm~^, appropriate for 
the outer regions of the disc. The density at a distance of a few grav- 
itational radii from the neutron star is expected to be > 10^^ cm^'', 
even for a radiation-pressure dominated disc ( Shakura & Sunvaev 
Il973 ^. which is beyond the range of validity for the reflection code. 
To illustrate the effects of higher density on the soft X-ray spec- 
trum, Figure 13 compares three reflection spectra computed with 
wh = 10^^ cm"'' with ones computed at 10^* cm~^, where the 
code remains valid. The figure shows very little difference in the 
spectra in the RXTE band including the Fe Kq line. This implies 
that the Fe Kq EWs computed from the tih = 10^^ cm~"^ models 
(Fig.|3} should still be useful guidelines for future observations. On 
the other hand, the three order of magnitude increase in density sig- 
nificantly alters the soft X-ray features predicted by the reflection 
spectra. The differences mainly lie at energies < 0.5 keV, where the 
continuum level is raised due to the increase in free-free emission 
and absorption. As a result, the EWs of many of the emission lines 
have been substantially decreased. The expectation is that increas- 
ing the density further will continue to raise the continuum level at 
soft energies, as more line emission is absorbed by the continuum, 
and decrease the EW of the remaining emission lines. Clearly, the 
soft X-ray emission features are a strong diagnostic for the accre- 
tion disc density. It would be therefore be very interesting to obtain 
rapid, time-resolved broad-band spectra of X-ray bursts to exploit 
the information in those spectral features. In the near future only 
the XRT onboard the Swift satellite will have these features, plus 
the necessary rapid response capability. Unfortunately, soft X-ray 
spectroscopy of burst sources will be challenging. LMXBs are old 
systems and many of them lie toward to the Galactic center, thus 
suffering from significant absorption due to neutral hydrogen along 
the line of sight. However, there are a number of LMXBs which re- 
side in globular clusters that may provide lower extinction columns 
and allow sensitive soft X-ray observations. Rapid follow-up of X- 
ray bursts from LMXBs globular clusters therefore provide the best 
chance of utilizing the full information obtained in the reflection 
spectra. 




Energy (keV) 

Figure 4. The solid lines show three reflection spectra computed with a 
density of nj] = 10^** cm""', while the dotted lines denote models with 
identical ^ but with nn = 10^^ cm^-^. A kT = 2.5 keV blackbody was 
the illuminating continuum for all cases. The increase in density results in 
only minor differences in the spectra at energies > 3 keV, but there are 
significant changes at lower energies. 
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